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Abstract 

By using an effective field theory for the electromagnetic interaction of spin waves, 
we show that, in certain antiferromagnets, the latter induce non-reciprocal effects in the 
microwave region, which should be observable in the second harmonic generation and 
produce gyrotropic birefringency. We calculate the various (non-linear) susceptibilities 
in terms of a few parameters the order of magnitude of which is under control. 
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1 Introduction 



The response of magnetic materials to electromagnetic fields gives rise to a reach variety of 
interesting phenomena In particular, non- reciprocal optical effects in antiferromagnets 
have received considerable attention during the last years 0, ^. The possible existence 
of certain phenomena, like second harmonic generation (SHG) or gyrotropic birefringency 
(GB), is dictated by the magnetic group of the given material, which may (or may not) allow 
suitable (non-linear) susceptibilities to be different from zero The optical wavelengths 

induce atomic transitions which provide a potential microscopic mechanism to obtain non- 
vanishing susceptibilities. Indeed, this is the case for the observed non-reciprocal effects in 
^^^203 0, 0. However, alternative mechanisms to produce such effects cannot be ruled out a 
priori, and may even become dominant at certain wavelengths. It is our aim to demonstrate 
that this is indeed the case for certain antiferromagnets when the electromagnetic fields are 
in the microwave region. This region is very sensible to collective magnetic effects which 
makes a field theoretical description appropriate. 

The low temperature low energy properties of antiferromagnets (with spontaneous stag- 
gered magnetization) are dominated by spin waves. The spin wave dynamics at low momenta 
and energy is very much constrained by group theoretical considerations 0. The symme- 
try breaking pattern SU{2) —>■ U{1) tells us that the spin waves must transform under a 
non-linear realisation of SU{2) [^. In addition to that the space group and time reversal 
must be respected by the dynamics. The continuum approach ensures that it is enough to 
consider the rotational part of the space group, namely the point group, and the primitive 
translations as well as time reversal. A systematic description of the spin wave dynamics 
fully exploiting the above group theoretical constraints has been provided in [§]. 

In this paper we apply the general framework described in ^ to work out the electro- 
magnetic response of certain antiferromagnets in the microwave region. We have chosen a 
crystal with no primitive translations mapping points with opposite magnetisations. The 
point group is taken to be 3m, but repeating the analysis for any other point group is 
straightforward. This choice is motivated by the Cr203 crystal, which shows interesting 
non-reciprocal effects in the optical region as mention before. We have calculated the linear 
and non-linear electric and magnetic susceptibilities, which turn out to depend non-trivially 
on the frequency of the incoming radiation. In particular non- reciprocal phenomena in the 
SHG as well as the GB are predicted to occur. These results apply to any antiferromagnet 
(with spontaneous staggered magnetisation) of arbitrary spin and crystal point group 3m 
such that the magnetic ions lie on the z-axis and no primitive translation mapping points 
with opposite magnetisation exists. 
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In order to simplify the notation we will take h = c = 1, which lead to a relativistic 
notation. So x = (t,x) and q = (ci;,k). Subindices fi = 0,1,2,3, where the first one 
represents the time component. Furthermore we work with holomorphic coordenates z = 
X + iy and z = x — iy. We distribute the paper as follows. In section 2 we briefly review 
some basic aspects of electromagnetic wave propagation in media in order to explain the 
appearance of non-reciprocal effects in SHG as well as GB. In section 3 we quickly review the 
framework described in full detail in p|. In section 4 we present the effective lagrangian. In 
section 5 we work out the effective action which describes the response to the electromagnetic 
field and give the (non-linear) electric and magnetic susceptibilities. Section 6 is devoted 
to a discussion. Finally, in the appendix we list all the terms which are not displayed 
sections 4 and 5 in order to make the presentation simpler. 



2 Electromagnetic Waves Propagation in Media 

It is the aim of this section to briefly review some features of electromagnetic wave prop- 
agation in media, which are relevant for the rest of the paper. In particular the phenomena 
of gyrotropic birefringence (GB) and second harmonic generation (SHG), in connection with 
non-reciprocal effects, which arise due to time reversal violation in the medium. 

Let us recall the Maxwell equations in insulating and chargeless media 

VD = 

V X H = ^oD 

V X E = -doB (2.1) 
VB = 0, 

which are to be supplemented with the constitutive equations 

D = E+P , , 

2.2 

U = B-M, 

where V and Ai are the electric and magnetic response of the medium respectively. V and 
A4 are functionals of the electric and magnetic fields and, of course, depend on the physical 
properties of the medium. From the two equations above we obtain 



V X V X E + 9o'E = - [d'oV + V X doM) , (2.3) 

which is going to be the basic equation in our discussion. Once V and A4 are given this 
equation describes the propagation of electromagnetic waves in the medium. 
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For definiteness, consider first the linear response of a liomogeneous medium to electric 
fields only. In this case, the most general form for the electric response is Q 



^{z) = J dxx^^z - x)E\x), (2.4) 

where x and z are space-time vectors. The tensor x^'^i^ ~ ^) depends not only on time, but 
also on space coordinates, which caracterises the spatial dispersive medium. We shall make 
the standard assumption that it varies slowly over the medium. This is equivalent to an 
expansion of the tensor x"'^{'^) k) in powers of k in momentum space. We obtain 

V\t, z)= J dt'x"\t - t')E\t\ z) + y dt'-^^^'it - t')d,E\t\ z) + ■ • • . (2.5) 

The first and second term in the electric response are the polarisation and the quadrupolar 
moment respectively. Notice that the tensors associated to each order of the multipole 
expansion depend only on the frequency of the electromagnetic wave. 

Let us see next how some terms in the multipole expansion of both the electric and 
magnetic responses give rise to qualitatively new observable effects. Consider (|2.5| ) together 
with the leading term in the multipole expansion of the magnetic response (magnetoelectric 
term) 

V = x'^'^E'' + 'y^''^dcE' 
M = a^^E 

Upon substituting these expressions in ( |2.3| ), the quadrupolar and magnetoelectric term give 
rise to the so-called gyrotropic birefringence as we shall show next. Consider the plane wave 
solution, E"(x) = E'^(g)e~*«^ + E"(-g)e^«^', of the equation (U). Then E^{q) fulfills 



(2-6) 



2 ™b 



E\q) = 0, (2.7) 



where = k°'/uj gives the propagation direction and its modulus the refraction index (recall 
that q = (ci;,k)). Suppose first that the quadrupolar moment, ■j"''"^, and the magnetoelectric 
term, a^"^, are zero. Non trivial solutions to this equation arise from the condition that the 
determinant of the matrix on which the electric field acts vanishes. The anisotropy of the 
permitivity tensor is generally responsible for this condition to yield two values of the 
refraction index for each propagation direction, which is known as birefringence. Namely, 
two different plane waves propagate in each direction with two different polarisations and 
two different velocities [0. If the quadrupolar and magnetoelectric terms are restored, they 
enter eq. ( |2.7|) through an effective permitivity tensor, e"'' — (^e°''^'^a'^'' — iuj^y"'^'^^ , with a li- 
near dependence on the propagation direction, which is known as gyrotropic birefringence 0. 
In particular, the equations governing the propagation in directions n and — n are different, 
which implies that the GB is a non-reciprocal effect, since these propagation directions are 
related by time reversal. 
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Let us consider next the non-linear response of the system to electric fields. Thus we 
have to add to ( p^) new quadratic, cubic, . . . , terms [Til] 



Viz) = J dxx^\z - x)E\x) + j dxdyx^^%z -x,z- y)E\x)E\y) + ■■■. (2.8) 

In particular, the quadratic term in the above equation leads to the appearance of second 
harmonic generation. Whenever exists in the field E°'{x) a contribution of frequency uj the 
electric response will have in addition two contributions, one of zero frequency and another 
one of frequency 2uj. Then, the electric response will be, in general, a superposition of plane 
waves with frequencies multiple of uj 

V\z) = VSiz) + V^Az) + VUz) + ■ ■ ■ , (2.9) 



leading, in turn, to a similar superposition for the electric field solution of eq. ( |2.3| ). The 
contribution to SHG comes from V2q{z) in the expression above, which to lowest order can 
be written as 

(2.10) 

P"(2g) = x''''%q,q)E''{q)E''{q). 

In order to be more explicit consider the non-linear expressions for the electric and 
magnetic responses, 

~ \ K (2-11) 

= iJ^^^^E^E^, 

Once they are introduced in (|2.3|) we have a complicated non-linear equation. However, the 
non-linear terms are usually small. Therefore, if we pass the linear part of the response to 
the l.h.s. we can calculate the electric field solution perturbatively, E = E(o) + E(i) + ■ ■ ■. 
E(o) is the solution of the homogeneous equation (eq. (p.7|) without the quadrupolar and 
magnetoelectric terms), i.e., a monochromatic plane wave of frequency uj. Then E(i) follows 
from the equation 



-2^ab _ ^ab ^ -a-bi ^b^^ (^)e-*?^ = -(2cu)' [x""" + Ti" e"^" fx''" \ E^q) (g)Efo) (g)e-*"''^ (2.12) 



abc I ^d^ade,,ebc 



It is clear that the solution of this equation requires in the l.h.s. an electric field of 
frequency uj = 2uj. This is called second harmonic generation. Notice, moreover, that the 
second term in the r.h.s. depends linearly on the direction of the wave number. Then when 
this term is non-vanishing we have non-reciprocal effects in the second harmonic generation. 

In the discussion above we have presented the simplest situations which lead to non- 
reciprocal effects. In magnetic materials, as the one we are interested in, V and M. depend 
both on the electric and magnetic fields. In this case, since B = n x E, the non-reciprocal 
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effects can be obtained from terms depending on tlie magnetic field in botli the electric and 



magnetic responses [jT2|, |T3|. Furthermore, the (generalised) susceptibilities are constrained 
by the magnetic points group of the crystal. Since we are considering GB and SHG, which 
are dynamical effects, only the elements without the time reversal operator of the magnetic 
group are to be considered ^. For Cr20s with the spins aligned in the third direction 
this is the 32 group. The allowed linear susceptibilities (relevant for the GB) are 

P' = XeE' + x^B' 

(2.13) 

M' = 7fE^ + 7f5" 

and the bilinear ones (relevant for the SHG) 

P^ = XeeE'E' + ^^ME^E' 

+ XebE'B' + xTbE'B' + xfsE^B^ 

+ XbbB'B' + 2x%'IB^B' 
p3 = 2xeeE"E~' + xfBiE"B~' - E^B') + 2xbbB"B^ 

(2.14) 

= YeeE'E' + 2YEiE'E^ 

+ '^'e'bE'B' + YebE'B'"^ + Yi^E^B' 

+ -i'i^B'B' + 2^%iB'B'' 
= 2^f^E'E~' + ^f^{E'B' - E'B') + 2-ff^B'B'. 

In the remaining sections we shall calculate the contributions to the generalised suscepti- 
bilities above due to the spin wave dynamics. We will start with a local effective lagrangian 
describing the interaction between spin waves and electromagnetic fields. Upon integrat- 
ing out the spin waves we obtain a non-local effective action for the electromagnetic fields. 



which is equivalent to having a free energy |12|, taking into account that Lint = —Hint- The 



electric and magnetic response, and hence all the (generalised) susceptibilities, can be easily 
obtained as follows 

pa^SSeff M^=^-^. (2.15) 

dE'' 'SB" ^ ' 
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3 Building Blocks 



In this section we present the basic building blocks in the construction of an effective 
lagrangian for the interaction between the spin waves and electromagnetic fields, and their 
transformations under the relevant symmetries. The method we follow was thoroughly des- 
cribed in a previous article |^. Here we shall only give a brief overview of it. 

As it was mentioned in the introduction the spin waves are the lowest lying excited states 
of the antiferromagnetic ground state associated to the spontaneous symmetry breaking 
SU{2) U{1). This tells us that the associated field, U{x), is an element of the coset space 
SU{2)/U{1) 10], which transforms under SU{2) as follows: 



U{x)^gUix)h^ig,U), 



(3.i: 



where g G SU{2) and h G U{1) is a local {U{x) dependent) element which restores gU{x) to 
the coset space. If the alignment direction of the local spin is the third direction U{x) can 
be written as 



U (x) = exp 



'-^[n\x)S' + n\x)S^' 

J -IT 



(3.2) 



where 7r'(x) are the spin wave fields. These fields in the complex representation have the 
form TT^ = {-K^ ± i7r^)/-\/2 and the generators are written as S± = ±. iS"^. 

In addition to the continuous SU{2) transformations the action must be invariant under 
the space-time transformations. In our case we take the Cr203 as the underlying crystal in 
which the spin waves propagate. Cr203 enjoys the crystallographic point group 3m. The 
transformation properties of the U{x) field under the 3m ® T elements are 





U{x) - 


g-iU{x)h\ 


I : 


U{x) - 


U{x)Ch\ 


ay : 


U{x) - 


-> g2U{x)h\ 


T : 


U{x) - 


U{x)Ch\ 



C = e 
Ct = -C. 



(3.3) 



The nontrivial transformation under / is due to the fact that this particular transformation 
maps points with opposite local magnetisation in the antiferromagnetic ground state. The 
primitive translations act trivially on U{x) and have not been displayed. 

The spin-orbit is an important interaction which produces a gap in the spectrum of the 
spin waves because it breaks explicitly the SU (2) symmetry. The breaking part is given by 
some additional terms in the Heisenberg hamiltonian 



h=y: J^JSSJ + E D.,(s. X s,) + E M^s^s: 



(3.4) 



<'j> 



<«j> 



<«j> 
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where the tensors D'^j and M^^ break the SU (2) symmetry. The order of magnitude of such 
tensors is D" ~ {Ag/g)J and M"'' ~ {Ag/gfJ, where for CrgOg, Ag ~ IQ-^^f |l5l. In 



order to introduce them in the effective theory we take their local limit and promote them 
to sources with proper transformations under SU{2). By combining these sources with the 
SU{2) generators we obtain objects which transform covariantly under SU(2), 

D,,^D;^S^ ^ gD,,g^ 
M = M''''{S^^S^ + S''(g)S^) {g®g)M{g^®g^). ^ ' ' 

Finally they must be fixed to their more general form compatible with the point group 
symmetry, namely, 



D,, = D-,S+ 



D- = -D± 

zz zz 



Dtz = -D^-z 



(3.6) 

M = M-+{S+ ® 5_ + S_ ® S+) + M33(53 ® ^3)^ 

Therefore the objects from which we construct our theory are the spin waves given by 
U{x), the derivatives, 9^, and the spin-orbit tensors, D^j and M^j'. Let us arrange them in a 
simple form which provides us elements with easier transformations properties under SU{2) 



U\x)idf,U{x) = a-{x)S+ + a+{x)S- + al{x)S^ 
U\x)D,,U{x) = d-{x)S+ + dUx)S. + dl(x)S' 



{u\x)®U\x))m(u{x)®U{x)) = m—{x){S+(g)S+) 

+ m++(a;)(5_ ® S_) 



+ m^%x){S^®S^) (3.7) 
+ m-+{x){S+0S^ + S^0S+) 
+ m^\x){S+^ S"^ + S"^ ^ S+) 



From ( p.l|) the transformation properties under SU (2) for the coefficients of the genera- 
tors are 

a-{x) e^^(^')a;(a;) 



o-uW al{x) + d^9{x) 
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m (x) — > e^*'"^^^m (x) 



dr,„(x) e''^^^^d„Jx) m ^(x) — > e^^^^^m 



dpg{x) dpg{x) m^^{x) m 



X 



(3.8b) 

{X) 

■m^^{x) —>■ ■m^^{x), 



i.e., the non-linear SU{2) transformation is implemented by a U{l)iocai transformation. The 
second transformation in ( pj.8| a) allows us to introduce a covariant derivative = 9^ ± ia^ 
acting on . Covariant derivatives acting on ds or ms are redundant and should not be 
considered (see [§|). 

The space-time transformations are given by 



a?, ^ a? 



e:{C3+,a,}: <| ^ (3.9a) 

^■■U}-- { d;^ ^ -rf^, (3.9b) 

^ab _^ ^ab 

where the symbols ^/i, and t/i represents the transformation of the subindex under the 
space and time transformations respectively together with the corresponding coefficient in 
each case; the d superindex is the complex conjugate of a. 

Next we present the way of introducing the coupling to the electromagnetic field. Since 
spin waves have no electric charge they couple to the electromagnetic field through the 
field strength tensor, i.e., direct couplings to the electric and magnetic fields. This kind 
of couplings does not break the SU{2) symmetry and in order to maintain the space-time 
symmetry we impose the field E transforms like a vector and the field B transforms like a 
pseudovector under 3m point group, whereas under time reversal these fields transform as: 

T : I ^ (3.10) 

Since the spin waves are fiuctuations of magnetic moments there exists another kind of 
coupling given by the Pauli term. The Pauli term breaks explicitly the SU{2) symmetry 
and a source with appropriate transformation properties must be constructed to implement 
its effect in the effective theory. In the Heisenberg lagrangian with the Pauli interaction. 
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written in the second quantisation language, 

L = ^ ^^{x,)i {do - z/iSB(x,)) ij{x,) + ■ ■ ■ , (3.11) 

i 

a source Aq{x) ~ /iSB(x) can be associated to the Pauh term, which transforms hke a 
connexion under time dependent SU (2) transformations, 

Ao{x) ^ g{t)Ao{x)g^{t) + ig{t)dog\t), (3.12) 

such that now the theory will be invariant under time dependent SU{2) transformations. 
Therefore the effect of the Pauli term is implemented in the effective theory by changing the 
time derivative by a covariant time derivative, 

do^ Do = do-tAoix), (3.13) 

and eventually setting Aq{x) = fiSB{x). Once this change is performed one has to keep in 
mind that the and aj] contain the magnetic field encoded in the covariant time derivative. 

At this point the two sources of electromagnetic coupling to spin waves have been con- 
sidered. 



4 Effective Lagrangian 

Now we are in a position to construct the spin wave interaction with the electromag- 
netic field for the antiferromagnet. The way we choose to do this is a perturbative one: 
the derivative expansion. Carrying out this expansion to a given order is meaningful for 
low energy and momentum with respect to the typical scales of the antiferromagnet, given 
respectively by the superexchange constant J ~ lOmeV and the inverse of the lattice param- 
eter 1/a ~ 0.1pA~^ (the velocity of propagation of the spin waves relates both parameters, 
V = J a. It has the following value in Cr20^: v ~ 10~^c |T3|). The characteristic energy 



and momentum of the system are given by the external inputs of the electromagnetic fields, 
which are the same, namely, u;, and therefore the space derivative is highly suppressed with 
respect to the time derivative, vdi ~ 10~^9o. The suppression of the spin-orbit tensor has 
already been given, D ~ 10~^J and M ~ 10~^J. Terms proportional to and M force 
the local magnetisation to be in the third direction an give rise to an energy gap ~ 10~^J 
for the spin waves. The amplitude of the electromagnetic field must be constrained for the 
expansion to make sense. First, we consider the Pauli term. Since it is associated to the 
time derivative it is suppressed by J. We will assume that the remaining couplings of the 
electromagnetic field come from vector and scalar potential minimal couplings in a micros- 
copic model. The former is associated to a link and hence suppressed by 1/ea whereas the 
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latter is associated to a time derivative and hence suppressed by J/e. Therefore the electric 
field will be suppressed by J/ea, whereas the magnetic field will be suppressed by 1/ea^. 
The microscopic model may also have non-minimal couplings to the electromagnetic field 
arising from the integration of higher scales of energy and momentum. These terms would 
be suppressed by the above mentioned higher scales and will be neglected. In any case, as far 
as they respect the SU (2) and crystal point group symmetries their only effect is to slightly 
modify the value of the constants in the effective lagrangian, which are anyway unknown. 

Any effect due to spin waves is expected to be enhanced when we approach their energy 
gap. This is why we shall choose the energy of the electromagnetic wave of that order of 
magnitude. When in addition the amplitude of the electromagnetic wave is tunned so that 
E Oq the following relative suppressions hold: 



vdi,eavB ~ IQ-^J. 

Once the above relation are given we are prepared to construct the relevant effective 
lagrangian, invariant under SU{2) and space-time transformations given by ( p.8| ), ( p.9| ) and 



The effective lagrangian at the lowest order in which electromagnetic field appears reads 



do, eaE, d 



m 



(4.1) 



( p.lOp , for the effect we want to study: Non reciprocal effects in SHG and GB. 





(4.2) 



-|- Zjm 



-+ 



+ Zsm' 



33 



+ ZA{dt,% - d-,-A)E' - {d-,a^ - dla,)E'] 



+ Ziot[{d+^% - d^^a+)E' - {d^-^a+ - d+.a^ )E'] 
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+ ZnE'E^ 

When we take into account the Pauh couphng in a^, contributions to SHG arise as 10^^^ 
and 10~^^ effects, and contributions to 10~^ appear in the case of GB. As it will be shown 
these contributions give rise to the desired non-reciprocal effects. Its is important to notice 
that if a primitive translation mapping points with opposite magnetisation existed the terms 
with a single time derivative above would not appear in the effective lagrangian. 

Our action has been constructed up to third order (10~^) and the contributions to non- 
reciprocal effects arise only from the Pauli coupling which is much more suppressed in (|4.1|) . 
Therefore we might expect other contributions at higher orders. This is indeed the case, 
but in order to keep manegable the number of terms in the main text, these remaining 
contributions are relegated to the appendix. 



5 Electromagnetic Field Effective Interaction 

Our purpose is to describe non-reciprocal effects in SHG and GB mediated by spin waves. 
Spin waves are responsible for an effective interaction of the electromagnetic field giving rise 
to susceptibility tensors where the properties of the material (spin waves) are encoded. 

Hence we realise that spin waves are not to be observed in this experiments and therefore 
they must be eliminated from our theory. The way to do this is by integrating them out in 
the functional (path) integral |jl6|] so that the new action depends only on the electromagnetic 
field. In order to perform the integration we have to write the action explicitly in terms of 
the spin waves. This is achieved by expanding (|3.2|), with the following result: 



S[n,E,B] 



dx 



TT TT 



doTT^doTT' 

U[doTT+{fiB' + XE') + don-{fiB' + XE')] 



1 



ifiUln+ifiB' + XE~') - n-{fiB' + XE')]B^ 



(5.1) 



+ ^fiXf^{E'B' + E'B' 



+ [hE'E' + hE^E'' 

where only the terms contributing to bilinear and trilinear electromagnetic fields in the 
effective action to be calculated are kept. These terms are the only ones needed to describe 
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the desired effects. The new constants which appear in ( |5.1| ) are combinations of those in the 
previous section. Although we do not know their precise values, their order of magnitude is 
fixed according to the counting rules given in (|4.1|). We list them below: 



hi ~ ieaf 



A ~ D 

, D 
A ~ ea— . 

u 



(5.2) 



Recall that D ~ 10"^ J stands for the size of the spin-orbit term. 

Notice that the contributions to non-reciprocal effects to the leading order (10~^^ and 
10~^^ for SHG and 10~^ for GB) come from terms with at most two spin waves, which 
permits us to perform a gaussian integration in the functional generator. In addition to this, 
it is worth mentioning that at the order given above the effects are produced at tree level, 
i.e., without loop contributions. 

Once the gaussian integration is carried out a perturbative expansion of the spin waves 
propagator in the presence of electromagnetic fields has to be made, considering the free spin 
waves propagator. 



a; e 



-iq{x-y) 



Piuj) 



(27r)4 ^ ' ' ' ' cj2-A2' 

as the unperturbed part, leading to the electromagnetic effective interaction lagrangian 



(5.3) 



S^ff[E,B] 



+ 
+ 



dxfl ^ 
dxdyfl 
dxdyfl 



hiE^'E' + 62^^ + -fxX{E'B' + E'B' 



-fiX [WdlP{x - y)B' + B^d^,P{x - y)E' 

- ^ifiX^E' (5^(9oP(x -y) + doP{x - y)B'^) E' 

1- 2x 

— -tn X 



(5.4) 



E' (^B'^doP{x -y) + doP{x - y)B^) B^ 
+ B' (B^doP{x -y) + doP{x - y)B^) E' 



The arguments of the electromagnetic fields have not been explicitly displayed. They must 
be understood as the nearest in the closest propagator. 

Given the transformations under time reversal ( p.lO| ) it is clear that the SHG, terms with 
three fields, presents non- reciprocal effects due to the interference of different terms. The 
same is true for the GB since bilinear terms proportional to the magnetic and electric field 
appear. 

From the action above, together with the additional terms given in ( [A.6|) , the electro- 
magnetic response of the Cr20^ due to spin waves, leading to non-reciprocal effects in SHG 
and GB, can be easily obtained using ( |2.13| )-( |2TT5| ). 
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The linear susceptibilities read 



Xf (a;) = 2bJ^ 

xfM = 2hf, (5.5) 

Recall that X^b{'^) ^"^^ 7lf (^) ^ise to the GB. Notice that here this effect is proportional 
to the gap of the spin wave spectrum, which is in turn due to the spin-orbit interaction. 
The susceptibilities contributing to SHG read 

Xf|(a;,u;) = Aei/,V[P(a;) - 4P(2a;)] + 2A/i/,^a;[PM - P(2u;)] 

+ hH,fyP{uj)[P{uj) + 2P(2c.)] 

Xfl(^, ^) = \xe2py[P{u^) + 4P(2u;)] - Xe^fy[P{uj) - 2P(2a;)] 

- \xhfMP{u) + 2P(2c.)] - h^d,flu'[P{u:)P{2u:)] 

- ^X^cflu^P{u)P{2u) - ?>hfluj 

xTsi^M = Xg,fMP{uj)-2P{2uj)] (5.6) 

Xfl(^, ^) = ^MV.'c^[P(c^) + 2P(2a;) - 6u;^P(a;)P(2a;)] 

- Xg,fluj\P{uj) + 2P(2cu)] - 6ji/> 
Xf|(u;,a;) = -2A(73/.'cuP(2cc;) + 2j2^'a; 

Xfl(^,^) = J//'A/2a;[P(^) + 2P(2^) - 6u;2p(^)p(2^)] 

Xf|(a;,a;) = -Xe^flu'P^u) + ]^Xe,flu^P{uj) 
Xf|(a;,u;) = -^A53/>P(c^) - 2j2/.'a; 

7|fi(u;,a;) = A^i/>P(u;) 

7|i(u;,a;) = -iA^?3/>P(c^)+j2/.'c. 

Yi^{uj,uj) = ^//2A/,2u;[P(a;) + 2P(2a;)-6u;2p(u;)P(2a;)]. 

Notice the non-trivial dependence in uo of the above susceptibilities. (This dependence 
is slightly more involved if one calculates the general susceptibilities x(^)t^') and ^({lo.lo') 
since the limit u = u' produces a few cancellations.) For this to be so it is crucial that 
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no primitive translation mapping points with opposite magnetisation exists. Otherwise only 
the local terms proportional to ji would survive. Notice also that the terms proportional to 
h arise due to the explicit spin-orbit breaking. Even thought this term gives rise to a local 
term in the electromagnetic fields at this order, it also contains explicit interactions with 
spin waves at higher orders. 

Let us finally mention that the terms proportional to ki in (|A.6| ) give rise to contributions 
to the quadrupolar momentum |P, |12| which are of the same order as the ones considered 
above. In fact, the quadrupolar terms give the unique contribution to SHG in a crystal 
which contains a center of symmetry as it is the case for Cr203 above the Neel temperature 
The associated susceptibilities can be easily calculated. They are local and will not be 



displayed exphcitly. 



6 Discussion 

We have used an effective field theory for spin waves in an antiferromagnetic material to 
describe its response to electromagnetic fields in the microwave region. The starting point 
is a local effective lagrangian which fully exploits the fact that spin waves are Goldstone 
modes of a SU (2) —>■ U (1) symmetry breaking pattern together with the crystal space group 
symmetry and time reversal. By integrating out the spin waves we obtain a non-local effec- 
tive action which encodes the response of the material to the electromagnetic field. From 
this effective action the various linear and non-linear electric and magnetic susceptibilities 
can be immediately obtained. We have given explicitly those relevant to the GB and SHG 
experiments. These susceptibilities depend on a relatively large number of unknown con- 
stants (~ 23) and a microscopic calculation is required to assign definite numbers to them. 
However, their order of magnitude can be readily established in terms of the typical lattice 
spacing a and the energy of the first gapped excitation J. Notice also that these susceptibil- 
ities present a rather non-trivial dependence on cj, the frequency of the incoming radiation. 
This dependence cannot be obtained from the magnetic group symmetries alone and it is a 
direct consequence of the existence of spin waves in an antiferromagnetic crystal where: (i) no 
primitive translation mapping points with opposite magnetisations exist, and (ii) spin-orbit 
effects are sizable. 

As mentioned in the introduction, from group theoretical arguments it has been known 
for long that certain antiferromagnets may support non-reciprocal effects However, 
group theoretical considerations alone are unable to indicate any mechanism by which these 
effects may be realised, not even to provide an order of magnitude estimate. In ref. |0] non- 
reciprocal effects were observed in Cr20^ in the optical region, and in ref. a theoretical 
explanation was presented. A microscopic mechanism leading to such effects was identified in 
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atomic transitions in which spin-orbit interactions and the trigonal field play a crucial role. 
Here, we have presented a totally different mechanism which also leads to non-reciprocal 
effects but in the microwave region. Namely, the interaction of spin waves with microwave 
radiation. Since the mechanisms are completely different, the susceptibilities given in |0 
hardly have any resemblance with ours, except for the general group theoretical constraints 
that both of them must fulfil. For instance, the susceptibilities in tipically depend on the 
size of the of the Cr ions, the energy differences between atomic levels, and matrix elements of 
perturbations (like spin-orbit terms) between atomic states. Instead, ours tipically depend 
on the lattice spacing, the Heisenberg coupling and the spin wave energy gap. However, 
some common features do exist. The spin wave energy gap is due to spin-orbit terms, which 
are crucial to obtain non- reciprocal effects both for ref. and for us. In fact, within our 
approach it is clear that spin-orbit terms are the responsible for the spins to point to the 
third direction, and hence for the magnetic group to be what it is. Also the magnetoelectric 
susceptibilities in ref. are proportional to the magnetic moment of the ion and so are ours. 
Nevertheless, for the optical region it seems that only the susceptiblities 7^;, xee and '-)ee 
are relevant and, then, it is crucial that xee ~ Iee 7^ for the observation of non-reciprocal 
effects in second harmonic generation. In our case all the susceptibilities (i.e. including xeb-, 
Xbb, Ieb and ■Jbb) have the same order of magnitude which provides further observational 
possibilities. For the sake of comparison, we give below our order of magnitude estimates for 
the susceptibilities given in assuming the orders of magnitude that we have been using 
for the parameters so far. From our •je, Xee and •jee we find, adopting the notation of 0, 
a^^ ~ 10~^ and x ~ 7 ~ 10~^^CN~^, which appear to be a few orders of magnitude larger. 

We have not been able to locate experimental results in the literature to test our formulas 
against. We expect them to become available at some point. It would be particularly 
interesting to be able to browse the microwave region with several frequencies so that the 
uj dependence in (|5.5| ) and ( ^.6|) could be checked and the free parameters fitted. If the 
incoming radiation is directed along the third axis then only 7|fJ(co', u), XeK^^j XsK^^j 
XEsi^y^) Xee^^^^) relevant. 

Although for definiteness we have focused on the Cr203 crystal, which has spin 3/2, our 
results hold for any antiferromagnetic crystal (with spontaneous staggered magnetisation) 
with crystal point group 3m and arbitrary spin, as long as no primitive translations mapping 
points with opposite magnetisation exist. This includes for instance V2O3 (spin 1). It is also 
worth emphasising that the method we have used is general enough to become applicable 
to any antiferromagnet of any spin and crystal point group, as long as there is spontaneous 
staggered magnetisation. The allowed terms in the effective lagrangian, however, depend on 
the particular crystal point group and on the particular distribution of the magnetic ions in 
the crystal. 
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Appendix A 

In this appendix we will present the higher order terms in the effective action (4.2) which 



contribute to he SHG and GB to the same order as those in (|0 
Contributions to fourth order (10"*^): 



i^a^DQaQ — DQaQ)E'^ (A. la) 



dlzidt-z^o + d^AW + dUd^A + rf+ flo 

dl-Adt-z^Q + d^-,al)E^ + dl,{d^,a^ + dt,a^)E~^ 

ML(4s«o + 4-f«o ) + dUd~z-A + d+,a^)]E^ (A.lb) 

ML(4.«o + d^A) + dUd-^A + 4i«o )]^' 
[dl-zAz^Q + c^3^4) + dl,{d^-,4 + dp,aQ)]E^. 

To fifth order (10^^°): 

(d+ Doao + d:^Do4)E^E^ + {dz,Do4 + d%D,a-,)E'E' 

{d+Doao + d^,Doa+)E'E' + {d^,D^a+ + dlDoa^)E'E' 

{dt^Doao + d-^Doa^)E^E^ + {d^^Doa^ + dtsDoao)E'E^ 

{dt,Do% + d^,Do4)E'E^ + (rfg- Do< + dlD^a^)E'E^ (A.2a) 

(c/+ ao + d;,a+)E'doE' + [d^.a^ + dt,a^)E^doE' 

{dt^a^ + d^^a^)E^doE^ + {djA + dt^a^)E'doE''^ 



^idtJ?, - d-^dt,){dlE-^E-^ - dl-X'E') 
^Azds, - d^Jt.){dlE^E^ - dl-X'E') 
lAA, - d:A,)idl-zE'E^ - dlE-'E-^) 
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z[(m+^d- - m-'dt,)E^E^ - {m-'dt, ■ 
^[{m+%, - m-'dt,)E^E^ - {m~'dl 



m+^dJ,)E'E' 
m+^d^^)E'E' 



^idtJ^,-d:Jt,)idlE-'E' 
t{dlds,-d^,dl){dlE'E' 
^(4.4-2 - d:Az)idl,E'E' 
i{dld~,,-d^,dt,){dl,E'E' 
i[(m+3rf- - m-^dt^)E'E^ - 
i[{m+^d^, - m~^dj,)E~'E^ - {m-^d+^ - m+^d^^)E'E^ 



dl.E^E') 
dl.E^E') 
dlE^E') 
dlE^E') 



[m 



-^dt 



m 



^'d-jE^E^ 



(A.2b) 



i{dl^E^B^ - dl,E-^B-^) 
i{dl,E^B^ - dlE'B') 
tidlE'B' - dl.E^B') 
i{dl-X'B' - dlE^B') 
iidlE^B-^ - dl.E'B^) 
i{dl,E'B'-dlE'B^). 

To sixth order (10~^^): 

^{dtJ~ - d;Jtz)iE"doE' - E'doE')E^ 

mtJs. - d;Jl)E'doE' - (dJ-Jt, - dt-J^,)E'doE^]E' 

iidtJs-. - d^^dt,){E^doE' - E'd,E^)E' 



(A.2c) 



(A.3a) 



{E'doE'' - E'doE')B^ 
E'doE^B'' - E^doE~'B' 



(A.3b) 



(E'd-,E' + E'd,E')E^ 
E'E'd^E^ 

E'E'{d,E' + dsE') (A.3c) 
E^E'\d,E' + d2E'). 

Where to reduce the number of terms in (|A.3| b) the homogeneous Maxwell equations, 
which are satisfied automatically, have been used. It is important to notice that if a primitive 
translation mapping points with opposite magnetisation existed, the terms in ( |A.l| a), ( [A.l| b), 
( |A.2| c) and ( |A.3| a) would not appear. 
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In spite of the large number of terms, we will see that most of them contribute in the 
same way to the effective action for the electromagnetic fields. 

Indeed, when we expand the terms, above in order to make explicit the interaction 
between the spin waves and the electromagnetic field, the terms below must be added to 
!T]) keeping bilinear and trilinear terms in the electromagnetic fields. 



AS[-K, E, B] 



dx 



+ idi{n+don+E' - n-don-E'') 
+ id2i7T+don-E^ - n-doTT+E'^) 
+ ieiU{doTT+E'doE' - do7r~E'doE') 
+ ie2U{do'K+dQE'E^ - doTT-doE'E^) 
+ tesMdon^E'doE' - don-E^doE^) 
+ ihU{7r+E'E'-n-E'E') 



TT~E'E^) 
TT'E'^B^) 

7i-E'B^) 
Tc'E'^B') 



(A.4) 



+ zf2M'^+E'E^ 
+ ig^U7^^E'B' 
+ ig2M7r+E'B' 
+ ig;U7^+E'B~^ 
+ ihfl{E'doE~' - E'doE')E'^ 
+ zjif^{E^doE'-E~'doE')B' 
+ ij2fl{E'doE^B' - E^doE'B') 
+ kifl{E'd-,E' + E'd,E')E'^ 
+ k2flE-'E~'d^E^ 
+ hflE^E\d,E^ + d,E') 
+ k^flE^E\d,E^ + d,E' 



Notice that no terms beyond two spin waves appear in the previous action. This does 
not change the procedure of gaussian integration carried out in section 5. The constants in 
have the following order of magnitude: 
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ea f D 

c ~ rfi ~ ea ( y 

e. - [-j) - (ea) 

gi ~ {eafDa h ~ {eafD"^ 

Ji ~ (ea)='^ /ci ~ {eafj. 



(A.5) 



The effective action for the electromagnetic fields in (|5.4|) must be augmented with the 
following terms 



ASeffiE, B] = I dxdyfil - -iXei [E'd'^Pix - y)E'E' + E'E'd^P{x - y)E' 



2 



1 



E'dlP{x - y){e2E'doE^ + e^doE'E^) 
+ {e2E'doE^ + e3doE'E^)diPix - y)E' 
'E'doPix - y)UiE'E' + fiE'^E"") 
+ {hE~'E~' + f2E'E^)doP{x - y)E' 
E'doPix - y){giE'B' + g2E'B^ + gaE^B'') 
+ {giE'B' + g2E'B^ + g^E^B')doP{x - y)E' 
j dxdudyfll^ - ^iX^E^{dlP{x - u){fiB^ + d2E''^)doP{u - y) 

+ doP{x - u){iiB^ + d2E'^)dlP{u - y))E' 
- ^i/i^A E' (dlP{x - u)B^doP{u - y) + doP{x - u)B^d^P{u - y)) B^ 

+ B' (dlP{x - u)B^doP(u -y)+ doP{x - u)B^dlP{u - y)) E' 
+ ^iX^cE^{dlP{x - u)E^dlP{u - y) + dlP{x - u)E^dlP{u - y)^E^ 



(A.6) 



+ -idiX^ 



E'[d^P{x - u)E'doP{u -y)- doP{x - u)E'd^P{u -y))E 



E'(dlP{x - u)E'doP{u -y)- doP{x - u)E^dlP{u - y))E' 
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dxf^ 



+ ijiiE^doE^ - E^doE^)B^ + miE^B'^ - E~'B')doE^ 

+ kiiE^d^E^ + E^d^E^)E^ + k2E^E^d^E^ 

+ k^E'E^id^E^ + d^E^) + k^E^E^id^E' + d^,E^) . 



These contributions to the electromagnetic effective action have been included in our 
final results for the generalised susceptibilities in formulas (|53|) and (|5.6|). 
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